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Abstract. We apply ideas from the cluster method to g-count the permuta- 
tions of a multiset according to the number of occurrences of certain generalized 
patterns, as defined by Babson and Steingrimsson. In particular, we consider 
those patterns with three letters and one internal dash, as well as permutation 
statistics composed of counting the number of occurrences of multisets of such 
patterns. Counting is done via recurrences which simplify in the case of per- 
mutations. A collection of Maple procedures implementing these recurrences 
accompanies the article. 



1. Introduction 

In its most basic form the cluster method, as described in 10J, counts the num- 
ber of words of length I in a given alphabet which avoid a certain set of forbidden 
subwords (i.e. contiguous blocks of letters). In that same article, Noonan and Zeil- 
berger extend the method to count the number of words of length t which contain 
a given number of those forbidden subwords. The permutation statistics des, inv, 
and maj were implicitly treated with this method in |13j . Here we consider distri- 
butions of generalized permutation patterns of length 3, as introduced by Babson 
and Steingrimsson [5j . Burstein and Mansour [B] lists the number of words of length 
I in the alphabet {1, .., k} avoiding a given generalized pattern for each pattern of 
length 3. Claesson and Mansour [Sj provide a recurrence to compute the number 
of permutations with r occurences of a generalized pattern of length 3 with one 
internal dash. In this paper, we derive recurrences to compute (via a generating 
function) the number of permutations of a multiset with r occurences of a general- 
ized pattern of length 3 with one internal dash. The special case of permutations 
of {1, ...,n} allows for some more specialized recurrences. The method can si- 
multaneously compute the multivariate distribution of different patterns, yielding 
recurrences to compute certain pattern-based permutation (multi)statistics. 

In section 2 we outline the conventions regarding operations on words and recall 
definitions of generalized patterns. Section 3 describes how we apply the cluster 
method to derive recurrences for the distributions of generalized patterns on mul- 
tiset permutations. Section 4 considers the special case of distributions of patterns 
over permutations in S n , adapting the recurrences found in Section 3. Section 5 de- 
tails the accompanying Maple program which provides several tools for the reader 
to verify and to explore the conclusions discussed in this paper. A summary of 
results and suggestions of future directions can be found in the final section. 
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2. Preliminaries 

2.1. Operations on Words and Notational Conventions. We will consider 
words w = W\W2 ■ ■ ■ we in the alphabet A = [n] = {1, 2, . . . , n}. Borrowing notation 
from the theory of formal languages, let A* denote the set of all (finite) words in 
the alphabet A. We will use t{w) or simply I to represent the length of word w and 
n to denote the alphabet size. In the case of permutations we use n to represent 
both. For words u — u\- ■ -u% and v — V\ ■ ■ ■ Vg,, we denote the concatenation as 
uv = u\ ■ ■ ■ UkV\ ■ ■ -V£. We use Vi(w) to denote the multiplicity of i in w, and in 

generating functions it will be useful to consider x w := J\ i>1 x v i ^ w \ where the Xj 
are commuting indeterminates. Note that x uv = x u x v . 

We may also define the reversal and complement operations on words in [n] . Let 
w r = wgwe-i ■ ■ ■ w\ and w c = (n + 1 — w\){n + 1 — w%) ■ ■ ■ (n + 1 — we). If we are 
given an alphabet vector m = (mi, . . . , m„) and consider the set of words such that 
Vi(w) — mi for each i, then reversal is a bijection. Complementation is a bijection 
if and only if mi = m„, m 2 = m n _i, . . . , m L „ /2J = mf n / 2 ] ■ 

Since we count the descents of a word, it will be useful to define the set Des(w) := 
{i : Wi > Wi+i}. This should not be confused with des(w) := \Des(w)\. 

2.2. Generalized Pattern Functions. Babson and Steingrimsson introduced gen- 
eralized permutation patterns in [2], giving definitions which apply equally well 
to words (i.e. permutations with repeated letters). For simplicity we work with 
words in the alphabet [n] — {1,2, ...,n}. A pattern is written in the alphabet 
{a, 6, c, . . .} where dashes may or may not separate successive letters. The ab- 
sence of a dash between two letters indicates that the corresponding letters in 
the word must be adjacent; otherwise corresponding letters may be arbitrarily far 
apart. Further the corresponding letters of the word must be order-isomorphic to 
the pattern. We will write the pattern as a function, where applying the func- 
tion counts the number of occurrences of the pattern. This is best understood 
through examples. The pattern function (a — bc)(w) for word w = w\W2 ■ ■ - we 
is the number of subsequences WiWjWj+i such that i < j and Wi < wj < Wj+\. 
Likewise (abc)(w) := \{i : Wi < lOj+i < w;+2}| and (ba) = \{i : Wi > tUj+i}|. For 
words we have additional (nontrivial) pattern functions with repeated letters such 
as (a — ab)(w) = ■ i < j,Wi — Wj < Wj+i}\. We will regularly use the 
symbols a and r to represent pattern functions. Treating these as functions allows 
multi-pattern functions like a + 2t, the number of occurences of a plus twice the 
number of occurences of r. 

In this paper, we will only explicitly consider those patterns with three letters 
and one internal dash between the first and second letters, that is those of the 
shape (x — yz). These are called type- (1,2) by Mansour [9]. If y > z, then we call 
such patterns descent-based and if y < z they are called rise-based. Patterns of 
type-(2,l) are easily handled via the reversal operation, since a r {w) = a(w r ) where 
(a: - yz) r = (zy - x). 

We are interested in the distribution of patterns over sets of words. For a pattern 
function a and set of words W, the distribution of a over W is the generating 
function 
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The W superscript is usually omitted. We will use ideas from the cluster method 
to derive recurrences which quickly compute F a over the set of all words in the 
alphabet [n], as well as the set of words with a prescribed number of each letter. 
Two statistics a, r are equidistributed over W if F^ (q) = F^(q). 



3. The Cluster Method for Recurrences 

We now adapt the cluster method to count the number of words with a given 
number of occurrences of a given pattern. Instead of counting the number of oc- 
curences of a descent-based [resp. rise-based] pattern, we instead consider the 
problem of weight-counting the descents [resp. rises], where the descent- weight 
[resp. rise-weight] is determined by the pattern. For example, the pattern (a — cb) 
can be weight-counted, 

(a - cb)(w) = : j < i,Wj < w l+1 < Wi}\ = ^ \{j : j < i,Wj < w i+1 }\. 

iEDes(w) 

Denote the descent-weight of i £ Des(w) according to pattern a by o~d(w, i). Hence 
(a - cb) d {w,i) = \{j : j < i,wj < tOi+i}|. 

Thus we may manipulate the sum in F a as follows: 



F a (q,x) = ^2 a <( 



= 2 xw n <f i{i ' w) 

w£A* i^Des(w) 

= E xW II (! + i^ w) - !)) 

w£A' SCDes(w)ieS 

where in the last step we expand the product. In essence, this is the principle of 
Inclusion- Exclusion (alternately, Mobius inversion or sieve methods). As in |10j we 
are counting pairs (w, S) which can be considered marked words, where some subset 
S of the descents of w are "marked" (as a teacher would mark a wrong answer). 
Each marked descent WiWi+\ has weight {(f^ w ' 1 ) — 1), and the weight of a marked 
word is the product of the weights of the marked descents. 

We break here to provide a concrete example. Consider the word w = 2637541 
and pattern a = (a — cb). Then ad(w,4) = 2, since 2 — 75 and 3 — 75 are both 
occurences of (a — cb) involving the descent at position 4. Similarly, <Jd(w, 2) = 1, 
o~d{w, 5) = 2, and o~ d {w, 6) = 0. Now consider the marked word (w, {2, 5}), better 
viewed as 2637541. Only the descents 63 and 54 contribute, yielding the weight 
(q 1 -l)(q 2 -I). Similarly, the weight of (to, {4, 5, 6}) is (q 2 - l)(q 2 - l)(q° - 1) and 
the weight of (id, 0) is 1 because of the empty product. It is easily checked that the 
sum of the weights each of the 2 4 marked versions of w equals q a ( w ) = q 5 

From the cluster method, we partition the set of marked words into three classes: 

(1) The empty word, which has weight q . 

(2) Those words ending with a letter which is not part of a marked descent. 

(3) Those words ending with a marked descent. 



4 



ANDREW M. BAXTER 



If a marked word ends with a marked descent, it must end in a marked descending 
run Wi > Wi+i > ■ ■ ■ > we, where each component descent is marked. We remove 
the maximal terminal marked run (also called a cluster), leaving a word in the 
second category above. To continue the example above, (2637541, (2, 5, 6}) has 541 
as its maximal terminal marked run. To contrast, note that (2637541,(2,4,6}) 
has only 41 as its maximal terminal marked run, since 54 is not marked. Since 
weights are multiplicative, the weight of (2637541,(2,5,6}) is the product of the 
weights of (2637541,(5,6}) (the cluster) and (2637,(2}) (the rest). Furthermore, 
rearranging the letters outside the cluster does not affect the weight of the cluster, 
e.g. (2637541, {5,6}) and (7362541, {5,6}) have equal weight. Thus we may fix a 
cluster and let the prefix vary, guiding our next transformations. 

For any subset T C [n], identify T = t\ > t 2 > ■ ■ ■ > tk with its corresponding 
descending run *i* 2 ■ ■ ■ tk, and so write wT — 10*1*2 •••*&■ Algebraically, the above 
partition translates into: 

F a { q ,x) = £ x w yi n^^- 1 ) 

w£A» SCDes(w)ieS 

= i + EE 1 ™ E n^^-i) 

i£[n]w£A* SCDes(w)ieS 

+ E E xWT E nfe"* 

TC[n],\T\>2w£A* SCDes(w) ieS 

/m-i 

= 1 + ]T xWT II (q ad(wTAw)+ ' 

TC.[n],TjtQweA* \ j=l 



Tl-l 



- 1) J| ( q <Td(wT,t(u>)+3) _ 1) 



3=1 



-i)| E Utf'^'O-i) 

SCDes(w) ieS 



To rephrase the preceding observation, ad(u>T,l(w) + j) is dependent only on 
the values of the Vi(w), and not on the order of the letters in w itself. For example, 
it is clear to see that (a — cb)d(wT,£(w) + j) = v\{w) + ^(w) + • • • + Vt j+1 -i(w). 
Therefore, let m = (mi, . . . , m n ) E N n and partition A* into classes A* m = {w G 
A* : Vi{w) — rrii for all i}. Note A* m is the set of multiset permutations with m 
giving the multiplicities of each element. Let dd{m, T, j) = <Td(wT, mi + • • • m n 
where w e A* m (it docs not matter which w you choose). This is the descent- 
weight of the j th descent in T in the word wT. For example, (a — cb) d (m, T,j) = 
mi + m 2 + ■ ■ ■ + mtj+i-i- Hence we can write F rT (q,x) as 

m-i 

F a (q,x)=i+ j2 ^■■■< n e xT n {cf d{m ' T ' j) -i) e e n^"^' - 1 ) 

roGN" TC[ti],T/8 j = l weA^ SCDes(w)ieS 

Let us focus on J2 weA , q 17 ^) (the coefficient of x™ 1 • • -x™ n in F a {q,x)), which 
we denote F a (m, q). We will write recurrences for F a (m, q) in the form of operators 
P a so that P a F a — F a . Let be the shift operator for nit, that is, 

E^fim-i, ...,m„.. .,m n ) = /(mi, . .. ,mj - 1, . . . ,m n ) 
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and Ej, 1 = ]J ieT E i 1 . Then it follows that 

T|-l 

F a (m,q)= E t X II (?^ (m ' T ' j) -l)^K?) 

TC[n],T#0 3=1 

for m — (mi,...,m n ) and each m, > 0. In the case that some m, = 0, then 
F a (m, q) = F a (m, q) for m = (mi, . . . , m,*_i, m i+1 , . . . , m„), since we may reduce 
by one each letter greater than i. Furthermore, if = for all i, then F a (m, q) = 1. 
These provide the initial conditions for the recurrence above. 

As we apply each E^ 1 operator, we will need to consider the shift's effect on 
crd(m,T, j). It is easily checked that 

E T 1 (a - cb) d (m,T,j) = mi + m 2 H h m tj+1 -i - \T\ + j + 1 

= (a^d)) d (m,T,j)-\T\+j + l 

since the descent tjtj + i follows all occurences of letters less than tj +1 except the 
\T\-(j + l) which lie in T itself. Similarly, 

E^ 1 (a - ba) d (m,T,j) = (a - ba) d (m,T, j) - 1 

E-\b^ba) d (m,T,j) = (b^ba) d (m,T,j) - 1. 

For a G {(ba), (b — ca), (c — ba)}, E^ x dd(rri,T ,j) — cfd(m,T, j). The operators for 
each descent-based a are summarized in the table below. 



Operators P a for descent-based a 



a 


Pa 


(ba) 
(a — cb) 
(b — ca) 
(c — ba) 
(a — ba) 
(b - ba) 


X/TC[n],2V0 ITj2l 1 (l 1 - l ) E T X 

V" 1 TTl T |- 1 /' " l i+m 2 H hmt i+1 -l—\T\+j+l 

LrC[n],T/a llj=l W J+ _ J-J-^T 

V -rTl T l _1 Cn" lt j+i+ 1+mt i+i+ 2H — nr 1 

2^TC[n],T#0llj=l W 3+ J+ J iJ-^T 

Etc [n] ,t/0 III, = i~ 1 (9™''' + 1 +mtj ' + 2+ ' ■ ■ +m " - 1)^1 

ETCM.r^nE'r 1 ^ 1 - 1 -!)^ 1 
ErcN^nPr 1 ^- 1 -!)^ 1 



For rise-based statistics, the process remains unchanged except now one considers 
marked rises (i.e. subsets of Rise(w) — {i : Wi < Wi+\}). To keep the relative 
ordering fo the ti G T consistent, we consider what happens upon the removal of 
the ascending run T r from wT r . 

F a (q,x) = 1 + £ Y, ^ ( II (<Z M ™ T " /W+3) -1) ) • 
TC[n],T^0 w£A* \ j=l J 

Following the same manipulations as for descent=based a, we obtain the corre- 
sponding operators shown in the table below. 

Since the reversal operator is a bijection in A* m , we see P a r = P a . Hence we 
have recurrences for F a (m,q) for each a of type-(2, 1) as well. Note that here the 
clusters are initial marked runs instead of terminal marked runs. 

Some operators appear multiple times among Tables [3] and [3] Since all patterns 
have the same initial conditions, we have proven the following equidistribution 
results. 
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Operators P a for rise-based a 



(ab) 
(a — be) 
(b — ac) 
(c — ab) 
(a — ab) 
(b - ab) 



Pa 



n/Jrv- w 

Erc[n],r^nJ= l r 1 («/ mi+ma+ - + ^ +i - 1 - 1)^ 1 

ErC[ti],T/( J|^l~ 1 (g" lt 3 + i+ 1 + mt j + i+ 2 + '" +mt i- 1 - 1)25-1 

Etc [„] ,7y n|2i~ 1 (g m S + 1 +mt ^ + 2 + ■ ■ - +m " 1 - l)^- 1 

ETcw^enSr 1 ^ 1-1 -!)^ 1 



Theorem 1. T/ie patterns {b — ac) and (b — ca) are equidistributed over A* m for 
any alphabet vector m. Similarly (a — ab) and (a — ba) are equidistributed over A* m , 
and so are (b — ab) and (b — ba). 

These are the only such non-trivial equi-distribution classes over multiset permu- 
tations among the patterns listed above, as evidenced by computing the distribution 
over A% 1 1 2 n (see table below). 



Distributions of descent-based statistics over A 



a 


F^q, (1,1,1,2)) 


(ba) 


1 + 18q + 33g 2 + 8g 3 


(a — cb) 


31 + 17g + llg 2 + g 3 


(b — ca) 


28 + 23g + 8g 2 + g 3 


(c — ba) 


37 + lOg + 9g 2 + 3g 3 + g 5 


(a — ba) 


60 


(b - ba) 


24 + 36g 


(ab) 


1 + 18g + 33g 2 + 8g 3 


(a — be) 


31 + 20g + 5g 2 + 4g 3 


(b — ac) 


28 + 23g + 8g 2 + g 3 


(c — ab) 


37 + 9g + 10g 2 + 3g 3 + g 4 


(a — ab) 


60 


(b - ab) 


24 + 36g 



(1,1,1,2) 



This method's strength lies in its applicability to multistatistics. Let a and r be 
two descent-based pattern functions and let Fr a ^(m,q,t) := J2 weA * q cr ( w )t T< - w \ 
Then just as before we may write F„ iT as 

F a , T (m,q,t)= ]T E HW^f'W-l) 

weA* m SCDes(w) i£S 

All of the above transformations hold, leading to the operator 

m-i 

P aT = J2 E T* II (q^ m ' TJ h f ^ m ' T ^ - 1). 

TC[n],T#0 j=l 

Of course this generalizes in the same way to multistatistics with any number 
of descent-based patterns of type-(l,2). Furthermore, specializing indeterminates 
(except x) allows one to consider compound statistics, such as the major index 
maj = (ba) + (a — cb) + (b — ca) + (c — ba) + (b — ba) + (a — ba) or multistatistic 
distributions like the Euler-Mahonian polynomials E?reS„ t maj( ^ q^^K Likewise 
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one may consider multistatistics involving rise-based pattern functions. The great- 
est restriction is that one may only combine descent-based statistics, or combine 
rise-based statistics, and must keep within type-(l,2) or type-(2,l). In short, one 
may combine patterns from only one of the four boxes below. 



(ba), (a — ba), (b — ba) 
(a — cb), (b — ca), (c — ba) 


(ab), (a — ab), (b — ab) 
(a — be), (b — ac), (c — ab) 


(ab), (ab — a), (ab — b) 
(be — a), (ac — b), (ab — c) 


(ba), (ba — a), (ba — a) 
(cb — a), (ca — b), (ba — c) 



Problems combining descent-based and rise-based patterns arise since both de- 
scents and rises are marked, resulting in non-monotone clusters. Thus for each 
T C [n] there is more than one cluster involving the letters in T . Furthermore one 
must consider submultisets T to compute distributions over multiset permutations. 
A variant of the above methods may still work for such multistatistics if one sacri- 
fices speed. As mentioned previously, the clusters for type-(l,2) patterns appear at 
the end of the word whereas clusters for type-(2,l) appear at the beginning. How 
to rectify this disparity eludes the author. 

4. Distributions over S n 

We now move to the special case of permutations, i.e. computing F a (m) for 
m = (1, 1, . . . , 1) G N n , which we will denote F a (n). In this case it is clear that 
the shift operator translates to the shift N~^ T ^F(n) = F(n — |T|). For each of 
(ba), (a — cb), (b — ca), and (c — ba) we will determine a recurrence to recursively 
compute 

F a (n) = <f W - 

These recurrences allow for polynomial-time computation, and the same method 
should allow for polynomial-time computation for multistatistics following the re- 
strictions discussed at the end of the previous section. In each subsection below, 
the cr subscript for F is omitted. Of course for each pattern a discussed below, the 
same recurrences will hold for the distributions of a r , a c , and a c r . 

4.1. The pattern (ba). From table [3] we get the recurrence 

F(n) = II (g'-Wn-m) 

TC[n],2V0 j=l 

Observe that qF(n) are the well-known Eulerian polynomials. 

4.2. The pattern (a-cb). From table [3] we get the recurrence 

n k — 1 

fc=l n>ti>--->t fc >l j—1 

since mi + m-i + • • • + rn tj+1 ~i = tj+i — 1. Define 

fc-i 

a(n,k)= Yl H(q ti+1+i - k -l) 

n>t 1 > - >t k >l j = l 
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so that F(n) = a ( n > k)F(n — fc), and let 

k 

b(n,k)= U 'l' ! * ' n - 

n>ti >--->t k >l j=l 

Conditioning on whether t% = n, we get that 

k — 1 k—1 

a(n,k) = £ n^ ,+1+J " fc - 1 )+ E IK ' ' ' 

n-l>ti>— >t fc >l j=l n-l>i 2 >-->tfc>l 3=1 

fc-1 

= o(n-i,fe)+ £ n^ 4 *' - *- 1 ) 

n-l>ti>— >is._i>l i=l 

= a(n — 1, fc) + b(n — 1, fc — 1) 
Similarly we can derive a recurrence for b(n, k) by conditioning on t\. 

k 

b(n,k) = II^"' 5 " 1 - 1 ) 

n>ti>->t fc >li=l 

e n(^ +j '- fc - i -i)+ ^ ( g «- fe -i)n(^+w-i) 

n-l>ti>--->t fc >l J=l n-l>t 2 >--->t fc >l i=2 

fe-1 

= 6(n-l,fc) + ( (Z "- fe -l) £ IK * 1: 

n—l>ti> •••>*&_! >1 j=l 

= &(n-l,fc) + (<? n - fe -l)&(n-l,fc-l) 

These recurrences, along with the initial conditions below yield fast computation 
of F(n). 

o(n,l)=n 6(n»l) = E2=i(9 i_1 -l) 

a(n, fc) = for k > n b[n, k) = for k > n 

This recurrence confirms the values listed in Table 2 in [5]. 



4.3. The pattern (b-ca). As in the case for (a — cb), we get a recurrence of the 
form F(n) — 5Zfc=i a ( n ' k)F(n — fc), where now we get 

fc-i 

a(n,fc) = £ J]^"*^ 1 " 1 - 1 )- 

n>tl>--->t fc >l 3=1 

We will also make use of a secondary function b(n, fc), defined as 

fc-i 

6(n,fc)= £ g-^n^'"^ 1 " 1 - 1 )- 

n>ti>->i fc >l j=l 
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Again conditioning on whether t\ = n, we obtain 

(fe-i 
£ ( g "-*i-l_l)JJ( g *i-*i + i-l_l) 
«-l>t 2 >— >tfc>l 3=1 

fe-1 fc-1 

= a^-l^ + g"" 1 X) g-^n^"^ 1 " 1 - 1 )- E J]^^ 1 " 1 

n— l>ta>-">t fc >l 3=2 n-l>t 2 >--->t k >l j=l 

= a(n - 1, ft) + q n ~ 1 b(n - 1, ft - 1) - a(n - 1, ft - 1) 

fe-l 

&(n,fc) = 6(n-l,fc) + ^ g^?" - * 3-1 - 1) JJfa 1 '-**- 1 - 1 - 1) 

n-l>* 2 >— >tjt>l 3=2 

(fe-1 
]T (g"-*-- 1 - 1) JJ^-*^- 1 - 1) 

n-l>t 2 >-->tfc>l 3 = 1 

= b(n - 1, ft) + 9"" (g n_1 6(n - 1, ft - 1) - a(n - 1, ft - 1)) 
With the initial conditions below we can quickly compute terms of F(n). 

a(n, l)=n b(n, 1) = EIU&T' - 1) 

a(n, ft) = for k > n b(n, k) — - for k > n 

Note that ParviainenpJ] has found closed-form formulas for the coefficient of 
q k in F(n) for 1 < ft < 8 and provides the automated method for higher k. This 
recurrence confirms the values listed in Table 3 in [5]. 

4.4. The pattern (c-ba). We again seek a recurrence of the form F(n) — Efe=i °( n ' k)F(n- 
fc), where 

fe-i 

a(n,fc)= J] n^* 3 - 1 )- 

n>*i>--->t fc >l 3 = 1 

We will also need to use the secondary functions 

fe 

b(n,k)= *fcll(9 n ~* , -l). 

n>ti>— >t k >l 3 = 1 

and 

fe 

n>ti>->t fc >li=l 

First observe that a(n, ft) = b(n — 1, ft — 1), since 

fc-i 

a(n,fc) = ^ ( tfc _i)JJ( g «-ti_i) 

n>tl>->t fc _l>2 3 = 1 

fe-1 

E ifell^ 1 ^- 1 ) 

n-l>*i>— >t fc _i>l 3 = 1 

= 6(n-l,ft-l) 



10 



ANDREW M. BAXTER 



Now consider the sum in b(n, fc) and condition on whether tk — 1- Then we get 

k k 

b(n,k) = y, -i)+ e n^'- 1 ) 

n>t!>--->t fc >2 j=l n>ti> — >t)s_i>2,t fc =l 3=1 

fe fe-1 

e (**+i)ii(« n - i - ti -i)+(9 n - i -i) e 1 

n-l>t 1 > - >t k >l j=l n— l>i 1 >»->i fc _i>l,t fc =l j=l 

= 6(n - 1, fc) + c(n - 1, fc) + (g"" 1 - l)c(n - 1, fe - 1) 
We also condition on to get a recurrence for c(n, fc). 

k fe — 1 

c(n,fc) = f[(l n ^- l )+ E (^ 1 -l)Il(^ tj - 1 ) 

n>ti> - >t fc >2 j = l n>ti>—>tk-i>2,tk=l j=l 

fc fe-1 

e n(^" i_ti - i )+(« n " 1 - 1 ) e n^"" 1 "*'- 1 ) 

n— l>ti>— >t k >l j=l n-l>ti>— >t fe _i>l j=l 

= c(n - 1, fc) + (g™ -1 - l)c(n - 1, fc - 1) 
This recurrence confirms the values listed in Table 1 in [S]. 



1) 



5. Accompanying Maple Package Cluster GPP 

A corresponding package of Maple procedures illustrating the above methods 
can be downloaded at this paper's website: 



http://math.rutgers.edu/"&axier/CZusierGPP/mde£./iimZ 

Load the file into Maple using the read command. The built-in help command, 
Help() will outline the program for you, describing the syntax for each of the main 
procedures (listed for your convenience in the table below. 

The reader is encouraged to experiment with the above package to see resulting 
distributions. 



Procedures in the Maple package ClusterGPP 



Procedure Name 


Description 


PatternCount 


Counts the number of occurences of a given pattern in a 
given word 


BFdist 


Computes the distribution of a given set of patterns over 
a given set of words. 


Rdist 


Computes the distribution F a (m, q) for a single pattern a 
using the operator P CT . 


RdistM 


Recursively computes the distribution of multiple patterns, 
e.g. Pcr iT (m) using the operator P a , T - This subsumes Rdist, 
as singletons may be entered. 



6. Conclusions and Future Directions 

The methods outlined above yield recurrences for the distribution of any number 
of multistatistics generated by certain combinations of generalized permutation 
patterns. Besides being interesting in their own right, setting the indeterminates 



APPLYING THE CLUSTER METHOD TO COUNT OCCURRENCES OF GENERALIZED PERMUTATION PATTERNS 



to in the distributions gives the corresponding pattern avoidance results, such as 
those results in |7;. For example, setting q = in the recurrence for the pattern 
(a-cb) above gives an alternate recurrence for the generation of the Bell numbers. 
Questions of pattern packing, as studied in [12l [IJ [4] , ask for the highest number 
of occurences of a given pattern or patterns. This is simply the degree of the 
distributions we have calculated. For example, the first 20 terms of the sequence 

a(n) := max{(a — cb)(ir) : ir 6 S n } 

can be computed in under a minute using procedures from ClusterGPP. The in- 
terested reader may adapt the above procedures to generate recurrences for the 
degrees of the polynomials. 

The methods above would need massaging to extend them to multistatistics in- 
volving both rise- and descent-based patterns, for example the multi-pattern func- 
tion (o — cb) + (a — be). One must overcome the hurdle of weight-counting both 
descents and rises simultaneously, and determining what to do with runs of marked 
descents and rises. To extend the methods to 4-patterns, such as (a — c — db), the 
weights on the descents would get much more complicated than they are above. Of 
course, certain combinations of 4-patterns are already amenable to our methods, 
such as (a — c — db) + (c — a — db) = (b — ca) + (a — cb). 
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